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The v — 12/5 fractional quantum Hall plateau observed in GaAs semiconductor wells is a suspect in the 
search for non-Abelian Fibonacci anyons. Using the infinite density matrix renormalization group, we find 
clear evidence that fillings v — 12/5 and v — 13/5 are in the k = 3 Read-Rezayi phase in the absence of 
particle-hole symmetry-breaking effects. The lowest energy charged excitation is identified as a non-Abelian 
Fibonacci anyon, distinguished from its Abelian counterpart by its local quadrupole moment. However, several 
experiments at v — 13/5 observe a re-entrant integer quantum Hall effect, implying particle-hole symmetry is 
broken. We rule out spin polarization as the origin of the asymmetry. Further, we point out extremely close 
energetic competition between the Read-Rezayi phase and a re-entrant integer quantum Hall phase. This com¬ 
petition suggests that even small particle-hole symmetry-breaking perturbations can explain the experimentally 
observed asymmetry between v — 12/5 and 13/5. We find that at v — 12/5 Landau level mixing favors the 
Read-Rezayi phase over the re-entrant phase. 

PACS numbers: 73.43.-f, 05.30.Pr 


I. INTRODUCTION 

The richness of the emergent excitations in many-body sys¬ 
tems can belie the simplicity of their underlying interactions. 
This precept underlies the continued effort to realize quantum 
materials that exhibit fractionalized non-Abelian quasiparti¬ 
cles. When a finite number of such non-Abelian quasiparticles 
are introduced to a system, they give rise to a set of degenerate 
energy levels which cannot be distinguished by any local ob¬ 
servable, and thus can encode quantum information resistant 
to decoherence. Such phases of matter, apart from their funda¬ 
mental interest, are the proposed building blocks of quantum 
computers resilient to decoherence . 1 2 

Two types of non-Abelian quasiparticles, Majorana zero 
modes and Fibonacci anyons, stand out for their potential ex¬ 
perimental relevance. While significant progress has been 
made towards realizing emergent Majorana zero modes in a 
variety of experimental systems,^ the number of candidate 
hosts for Fibonacci anyons—which are in some sense an in¬ 
teracting generalization of Majoranas—remains much more 
limited. When two Fibonacci anyons V approach each other, 
either one or both of them is annihilated: r x r -A 1 + r. 
These two possibilities can be used to encode a qubit of infor¬ 
mation, which is naturally protected from decoherence when 
the pair of r are spatially far separated . 1 The Hilbert space 
size grows with the number of Fibonacci anyons, following 
the Fibonacci sequence. Unlike Majoranas, Fibonacci anyons 
have “universal” braiding statistics: braiding Fibonaccis alone 
is sufficient to approximate any quantum gate acting on their 
space of degenerate states . 9 Various lattice models have been 
proposed realizing Fibonacci anyons, but all such models re¬ 
quire complex many-body interactions! HKI 

Remarkably, there has long been a suspicion that Fibonac¬ 
cis already exist as the low-energy excitations of a fractional 
quantum Hall state in the v — ^ plateau of GaAs quantum 
wells . 15 18 Experiments by Kumar et al^ have observed an 
incompressible state at v = ^ with a gap of about 80 mK, 


larger than expected from the model of non-interacting com¬ 
posite fermionsP 9 A theoretical proposal for a novel phase 
at v = yy and ^ filling was first put forward by Read and 
Rezayi.^ They described a class of incompressible phases— 
the Zfc “Read-Rezayi” sequence—which occur at filling frac¬ 
tions v = N ± ^qp2- The k = 3 member (“Z 3 state”) at 
filling v = 2 + | and v = 3 — |, which we call RR 3 and 
RR 3 respectively, involve “pairing” of triplets of particles and 
support Fibonacci excitations. 

There are potentially less exotic explanations of the v = ^ 
plateau }® 23 and one worries nature may operate with a prin¬ 
ciple of parsimony. While a number of interferometry ex¬ 
periments have provided suggestive evidence for non-Abelian 
braiding,®^ direct unambiguous detection of these anyons 
remains challenging. Thus realistic numerical simulations 
play an important role. Refs. [2 0 and [27lpresented some ev¬ 
idence that the ground state of the Coulomb interaction at 
v =s ^^ is described by the RR 3 phase. These exact diag- 
onalization studies were limited to small systems and required 
two assumptions: ( 1 ) the electron spin is fully polarized; and 
(2) the cyclotron energy is infinite (i.e., the absence of ‘Lan¬ 
dau level mixing’). Taken together, these assumptions cannot 
account for a striking experimental observation. In contrast 
to the v = ^ plateau, the v = ^ plateau is replaced by a 
reentrant integer quantum Hall state (RIQH) with quantized 

2 

Hall conductance a xy = 3 ^, which is believed to imply the 
formation of charge density order (CDO)PSESE3[lSl y e t under 
assumptions ( 1 ) and ( 2 ), fillings ^ and ^ are related by a 
particle-hole symmetry. Thus a compelling numerical case 
for Fibonacci anyons at must drop these assumptions in 
order to account for the RIQH phase observed at v = A/. 
Furthermore, it has thus far been impossible to measure the 
energy and nature of the anyonic excitations. These are the 
main aims of this work. 

In this work we use recent advances in our understand¬ 
ing of quantum entanglement and the infinite density matrix 
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renormalization group (iDMRG) to determine the nature of 
the v = ^ and v = ^ plateaus. These developments al¬ 
low us to measure the properties of individual anyonic exci¬ 
tations and to relax assumptions ( 1 ) and ( 2 ), and our findings 
are as follows. After verifying the existence of the RR 3 phase 
under assumptions (1) and (2) (Sec. 0 , we compute for the 
first time the energies and charge distributions of the anyon 
excitations (Sec. III). We find that the lowest energy charged 
excitation is a Fibonacci anyon, with a (disorder free) charge 
gap of approximately A ^ 2K. Fortuitously, the two Fi¬ 
bonacci fusion outcomes have strikingly different charge dis¬ 
tributions, implying the local quadrupole moment could be 
used for read-out of the fusion outcome. We then relax the 
assumption of spin polarization ( 1 ) by simulating both spin 
species, and verify that the system spontaneously polarizes at 
both v — ^^ in the absence of external Zeeman field. Fi¬ 
nally, we find that there is exceptionally close competition be¬ 
tween the RR 3 phase and a RIQH phase; in fact, changing the 
width of the quantum well drives a phase transition between 
the two, which may explain the absence of a v — ^ plateau in 
narrow wells. Previous studies likely did not detect this tran¬ 
sition because the CDO of the RIQH is strongly frustrated on 
the traditionally-used sphere geometry; exploring this com¬ 
petition is well-suited to the infinitely-long cylinder geometry 
used in iDMRG. By relaxing assumption (2) and including the 
effect of Landau level mixing, we show that the RR 3 becomes 
preferred over RIQH at v = ^, while within the limits of our 
numerical accuracy at v = it does not. This provides a 
plausible explanation for the experimentally observed asym¬ 
metry, further strengthening the interpretation of the v = ^ 
as the RR 3 phase, and sheds light on the probable nature of 
the v = ^ RIQH phase. 


II. IDENTIFICATION OF v = ^ AND f QUANTUM 
HALL STATES AS THE k = 3 READ-REZAYI PHASE 


We study electrons in a strong magnetic field Bz interacting 
via the Coulomb potential: 31 


T^Couiomb _ [ e 2 :p(r)p(r'): 

J r,r / 47T6 |r - r'| 


( 1 ) 


where r = (x,y,z) is a position vector. The electron gas is 
confined in the ^-direction by an infinite square well potential 
of width w, and so we project into the lowest subband of the 
square well. 3132 In a magnetic field the single particles states 
organize into Landau levels (LL) separated by cyclotron en¬ 
ergy huo c . In this Section, we assume the Coulomb interaction 
is projected into the spin-polarized N = 1 Landau level (LL). 
We will later enlarge the Hilbert space by including both spin- 
species and additional LLs. Energies are expressed in units of 
E c = e 2 /47r els, where £b is the magnetic length. For ref¬ 
erence, a 5 T magnetic field corresponds to Eq ~ 120 K for 
GaAs samples. Without internal degrees of freedom (spins 
or LL indices), the ^ and ^ states are related by an exact 
particle-hole symmetry, so to begin we only present results 
for v — ip 



FIG. 1. Orbital entanglement spectra for the two degenerate ground 
states at circumference L — 32 £b, well-width w — 3 £b- For each 
ground state, we pick two orbital cuts and plot the entanglement en¬ 
ergies vs. angular momentum. The low-lying spectra (highlighted in 
red) agree with the CFT prediction given in Eq. (2j), and thus provide 
an unambiguous identification of the RR 3 phase. 


We use infinite DMRcPSl to study the Coulomb interac¬ 
tion 0 on an infinitely long cylinder of circumference L, 
where its ground state is expected to have topological degen¬ 
eracy. Rather than trying to identify the phase via the ground 
states’ overlap with trial wavefunctions, we examine patterns 
in the quantum entanglement of the ground states which serve 
as the “order parameter” for topological order. For the RR 3 
phase, we expect ten degenerate ground statef 20 which split 
into two groups of five; within each quintuplet the states are 
related by translating the center of mass of the particles. We 
label the representative ground states from each quintuplet as 
|Qi) and |Q r )P^ 

The first evidence for the RR 3 phase comes from the en¬ 
tanglement spectra (Fig. |T|). We partition the infinite cylinder 
into left (L) and right (R) halves, each semi-infinitely long. 
Given a wavefunction |\h) on the cylinder, the entanglement 
spectrum 36 {e a } is the set of eigenvalues of — log pL , where 
Pl = Tr r 1 4/) (4/ 1 is the reduced density matrix for the left 
half of the system. Each eigenvalue e a is called an “entangle¬ 
ment energy” level, and carries quantum numbers of charge 
and angular momentum, just as a physical edge of the cylin¬ 
der would. Generically the low-lying levels of the entangle¬ 
ment spectra (along with their quantum numbers) mimic the 
physical edge theory of the phase,which can be used to 
identify the topological order. For the RR 3 phase, the cor¬ 
responding edge structure is a product of the Z 3 parafermion 
conformal field theory (CFT ) 39 with a U(l) boson. Each en- 
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FIG. 2. (a) The splitting of the topological degeneracy i£(|0 T )) — 
i^dQi)) as a function of cylinder circumference L is consistent with 
an exponential decrease, w is the width of the well, (b) Scaling of the 
momentum polarization M with circumference squared L 2 , which 
reveals the shift 8. Dashed lines indicate the theoretical predictions 
for the 1 and r ground states of RR3, and dotted lines the competing 
hierarchy (AH) and Bonder son-Slingerland (BS) states. 


tanglement spectra corresponds to one of six primary fields of 
the Z 3 parafermion CFT: 1, 'ip, ip^, e, a, aK The predicted 
level countings (i.e., the number of low-energy states for each 
momentum) are 


1: 1,1,3,6,12,..., 
V’M = 1,2,5,9,18,..., 


e : 1,3,6,13,24,..., 
(j/(7 ] : 1,2,5,10,20,.... 


( 2 ) 


Figure[l]shows the orbital entanglement spectra for the ground 
states |fii/ T ) at L = 32 £b- The pattern of the low-lying levels 
(highlighted in the figure) is indicative of a chiral edge mode, 
with level counting consistent with the corresponding CFT of 
the Read-Rezayi phase. In each spectra, the first four levels 
unambiguously match the theoretical prediction. 

Further evidence for the RR 3 state comes from the vanish¬ 
ing splitting of topological degeneracy in the thermodynamic 
limit, as shown in Fig. |2(a) For circumferences L > 17 £b> 
we consistently find ten ground states, as expected for the 
RR 3 phase.® Moreover, the splitting of the degeneracy is con¬ 
sistent with an exponential decrease with system size. 

Finally, a robust quantitative evidence for the RR 3 phase is 
given by the “momentum polarization,’E 1 43 which effectively 
computes the modular T-transformation (Fig. |2(b)| ). The mo¬ 


mentum polarization M is defined to be M = Tt[plK] with 
K the angular momentum operator on the cylinder; it mea¬ 
sures the average amount of momentum carried in the left half 
of the system. As explained in Ref. [Ml M encodes three topo¬ 
logical invariants: the “shift’® 8 , the chiral central charge c, 


xIO ' 3 



FIG. 3. The density profile of the charge-e/5 anyons in the pres¬ 
ence of a Gaussian pinning potential of width « 4 £b, computed 
on a L = 21£b cylinder. We enforce the presence of an any on 
in the central region by using various ‘topological boundary condi¬ 
tions’ [l, r], and minimize the energy subject to this constraint, (top) 
The Abelian g any on. (middle) The Fibonnaci § r any on. (bottom) 
When using topological boundary conditions [r, g r], the fusion rule 
t x r = 1 + r implies that there are two distinct anyon types that 
could arise: either an Abelian g or the non-Abelian g r. We find the 
non-Abelian type has lower energy. Consquently, the local properties 
of the excitation (e.g. the charge density and energy) agree between 
the (middle) and (bottom). These quasiparticles are about 15 Ib in 
diameter. 


and the anyon topological spin h a : 




vSI? 

( 4 n £ B ) 2 + 


^+0(e-i). (3) 



with the RR 3 phase and definitively ruling out the v = | 
hierarchy phase (8 = 6 ) and the Bonder son-Slingerland (BS) 
phase (8 = 4). 2 ^® The intercepts h a — ^ are also consistent 
with RR 3 (cf. Appendix [a|. 


III. THE NATURE AND ENERGETICS OF THE ANYONIC 
EXCITATIONS 

Having established the properties of the ground states, we 
turn to the anyonic excitations. While certain properties can 
be inferred from theory—e.g., the anyon types, charges, and 
fusion rules—we are interested in non-topological aspects 
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which depend on microscopic details. We first determine 
which anyon type is the lowest energy charged excitation. The 
gap to this excitation controls the stability of the phase, and is 
measured in activated transport measurements. Second, we 
measure the charge density profile of the anyons; if distinct, 
the two fusion outcomes r x r = 1 + r could potentially be 
detected by an electrostatic measurement. 

There are ten anyon types in the RR 3 phase. For each 
fractional charge m| (m G Z 5 ), there is one Abelian and 
one Fibonacci type, which we label np and respectively. 
The Coulomb interaction will generally drive charges to frac¬ 
tionalize into the smallest possible units, with charge ±|, but 
there remains two possibilities: ±| and d=|r. If the lowest 
energy excitation is =b| r, then charged defects would natu¬ 
rally nucleate Fibonacci anyons; braiding and fusion of the 
| r excitations would in principle be sufficient for universal 
topological quantum computing. 

We use the recently developed “defect-DMRG” method® 
in order to trap a single anyonic excitation 4 a’ on an infinite 
cylinder. Referring to App.[F]for details, we take advantage 
of the DMRG variational ansatz to ensure that to the left and 
right of some large central region, the state relaxes to one of 
the degenerate ground states l and r respectively; we call the 
lowest energy state with this ‘topological boundary condition’ 
|[Z,r]). The anyon a trapped in the central region must ap¬ 
pear in the fusion product: r x l = a + ..., where l is the 
anti-particle of l. In particular, when 1 = 1 = 1 is the bound¬ 
ary condition corresponding to the trivial anyon, the trapped 
anyon a is fixed by the right boundary condition, a — r. By 
measuring the energy of this configuration relative to the vac¬ 
uum, we determine the energy of the anyon a. The experiment 
can be done both with and without an electrostatic trapping 
potential. 

Beginning without a trapping potential, we define the en¬ 
ergy E a of anyon a to be the energy required to add a to 
the groundstate in the absence of a pinning potential plus 
the electrostatic interaction between a charge Q a point-charge 
and a neutralizing medium. With this definition, we find 
Ee = -0.0508 and E %r = -0.0511 at L = 21£ B ,w = U B . 
The energy of the Fibonacci particle is lower for all well- 
widths w = 0-3 ^ B and for both charges dbf. We show this 
energetic difference (at L = 21 i B ) in Fig. [4] Thus, in the 
absence of a pinning potential, it is energetically favorable to 
bind a Fibonacci r to a d=| charge. 

The “charge gap” for the RR 3 phase is the energy required 
to create and separate a pair of charge + §, — § particles from 
the ground state, A = E+ e T + 5_ g r , which we find numeri¬ 
cally to be about 0.017. Note that while the chemical potential 
/jl (which we have set to zero) factors into each E± g r individ¬ 
ually, the dependence cancels in A. This corresponds to about 
2 K at 5 T, much larger than the 80 mK activated-transport gap 
observed in experiments . 16 A similar discrepancy was found 
for the Moore-Read state at v = |, where numerics find a 
charge gap of about 0.022 (« 2.5 K at 5 Tj® while the ex¬ 
periments of Dean et al. 47 report an activated-transport gap of 
about 0.0047 (« 540 mK at 5T). At v = § the discrepancy 
is reduced somewhat by including the effects of Landau-level 
mixing,® but it is believed that disorder broadening plays a 


large role as well, for example an estimate of T & 1.75 K in 
the sample of Kumar et al. 16 Furthermore, we will show be¬ 
low that the size of the charged excitations is comparable to 
the length scale of the disorder, which could complicate the re¬ 
lation between the our clean charge gap and the gap observed 
from activated transport. 

To localize the anyons, we include a weak Gaussian pinning 
potential of width 41 B . The charge density of the | and | r 
anyons is shown in Fig. [3] Their charge distributions are qual¬ 
itatively distinct, so in principle they could be distinguished 
by their quadrupole moment Q zz , where £ points normal to 
the plane. Because of the concentrated charge distribution of 
the Fibonacci anyon, it receives a significantly lower poten¬ 
tial energy from the pin than the Abelian charge. This sug¬ 
gests an array of Fibonacci anyons can be pinned by trap¬ 
ping ± |-charged particles with a weak one-body potential 
and slowly cooling the system. Such a setup may be employed 
in a measurement-only quantum computing scheme.® 


IV. EVIDENCE FOR SPONTANEOUS 
SPIN-POLARIZATION 

If the system does not spin-polarize at either v = ^ 
or the two plateaus are not related by any symmetry, 
since the spinful particle-hole symmetry exchanges 2 + z> o 
4 — v. There have been few experimental studies on the spin- 
polarization of the ^ plateau. 49 The experiments of Ref. [49] 
found that applying an in-plane 5-field, which increases the 
Zeeman splitting, drove the ^ plateau through a transition; at 
a critical in-plane 5-field the gap measured from activated 
transport closed, then reopened at larger field. One expla¬ 
nation, believed to explain similar physics at v = || EE] 
is that the Coulomb point is spin-unpolarized until a criti¬ 
cal effective Zeeman field spin-polarizes the phase. Ref. 49 
pointed out that the spin transition observed at v = ^ re¬ 
sembles to some extent the transition at v = |, which likely 
occurs between the 332-Halperin state 21 and the Abelian hi¬ 
erarchy (composite fermion) state. However, the transport is 
anisotropic after the transition, and tilt-field experiments are 
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FIG. 4. Anyon excitation energies Eq , a at various quantum well 
thicknesses, for the RR3 phase in absense of disorder. (Left) The 
difference between the Abelian and Fibonacci anyons for charges 
Q = =b |. The data shows that for a fixed charge Q, it is energetically 
favorable to trap a Fibonacci rather than an Abelian anyon. (Right) 
The energy require to disassociate a pair of =b| quasiparticles. See 
the main text for caveats in regards to this data. 
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complicated to interpret because the in-plane B field com¬ 
bines with the finite well width to change the interaction in 
important ways, for example by mixing in higher sub-bands 
with a N = 0 LL character and reducing the effective well- 
widthP 

Rather than attempting to model the experiment in Ref.l49l 
we determine whether the ground state of the Coulomb Hamil¬ 
tonian 0 spontaneously spin polarizes when the Zeeman 
splitting vanishes. Indeed, the typical Zeeman splitting is 
small in comparison to the Coulomb energy, q^bB/Eq ~ ^, 
hence the spin-singlet vs. spin-polarized character of a plateau 
is largely determined by interactions. 

We check the spontaneous spin-polarization of the 
Coulomb state by keeping the full Hilbert space of both spin 
species with valence density v — jr + | and enforcing num¬ 
ber conservation of each spin separately (we ignore fully filled 
v = 2 and LL-mixing). At w = 3 £b and L = 21 £b we ob¬ 
serve long-range ordering of the fermion spin in the XY plane, 
signaling spontaneous breaking of SO(3) (see Appendix [c]). 
Despite the larger Hilbert space, the energy obtained agrees 
with the spin-polarized filling to excellent precision. For small 
system sizes that can be studied by exact diagonalization we 
also find that the ground state as well as the low-lying energy 
spectrum are fully spin polarized. In contrast, an identical nu¬ 
merical experiment at v = | is clearly spin-unpolarized, with 
an energy difference of 0.002 per flux below to the polarized 
phase, in agreement with experiment. 

It would be strongly desirable to revisit the experimental 
issue with a probe other than tilt-field, such as Knight shift, as 
proved effective at v = | ESE31 


V. THE RIQH PLATEAUS AND THE f, f ASYMMETRY 

In several experiments, the filling v = ^ lies within a re¬ 
entrant integer quantum Hall plateau (called “R2c”) with Hall 

conductance a xy = 3 ^. 28 29 Since Galiean invariance implies 

2 

a xy = translation invariance must be strongly broken. It 
is thought that charge-density order spontaneously develops 
and is pinned by disorder, rendering inert the fractional filling 
of electrons/holes in the valence LL. A similar RIQH plateau 
(“R2b”) with <j xy = 2^- lies directly proximate to the v = ^ 
plateau. Ref. [29] found that the partial filling at the center of 
the RIQH plateaus, v^b, ^ 2 c , very nearly satisfies particle-hole 
symmetry: 1 — c = ^ 2 b — 0.006. However, the width of the 
R2b plateau is a bit thinner than R2c, so v — ^ filling just 
escapes the RIQH region and develops as a separate plateau.®! 
Thus even small particle-hole breaking effects could account 
for the ^ ^ asymmetry by shifting the RIQH-RR 3 phase 
boundary. 

The nature of the charge density order in the RIQH R2b 
and R2c phases is not known a Hartree-Fock (HF) 

analysis predicts that among mean-field states either a two- 
particle “bubble” (a Wigner crystal of two-electron droplets) 
or a stripe (smectic) phase may be competitive at 2/5 and 
3/5 partial filling. 57 ®2l For a pure N = 1 Coulomb interac¬ 
tion the stripe is predicted to occur at a wavelength Ahf ~ 


4.91^b® It is expected that effects beyond HF will spon¬ 
taneously modulate the density along the stripes, effectively 
forming an anisotropic bubble phase without the symmetries 
of the triangular-lattice 


A. Charge density order 

To explain the asymmetry between v — ^ and , it is nec¬ 
essary to find the competing charge density order of the RIQH. 
CDO is subtle to study in finite-size numerics: the sphere ge¬ 
ometry will strongly frustrate CDO, while the torus has to be 
tuned to the correct aspect ratio (this is further discussed in 
Appendix [E]). On the infinite cylinder we must consider the 
finite circumference L and the unit cell used in the infinite 
DMRG. At filling p/q the unit cell must contain m • q flux 

2^2 

(m £ Z), corresponding to a period A m = m • q —jf- along 
the length of the cylinder. The RR 3 phase is commensurate 
with m = 1 for all L. 

In the numerics presented so far, CDO was implicitly for¬ 
bidden because we used the minimal DMRG unit cell m = 
1. To remedy this, we run simulations (still neglecting LL- 
mixing) for a variety of circumferences L and unit cells m • q 
chosen to accommodate the bubble and stripe phases, even in¬ 
cluding pinning potentials to preference various orders in the 
initial stages of the DMRG. For Coulomb interaction and fi¬ 
nite well width, the only CDO found is a stripe phase with a 
wavelength Acdo quite close to Ahf- Note that even though 
our geometry is quasi-ID, spontaneously breaking translation 
is not forbidden by Mermin-Wagner because the magnetic al¬ 
gebra renders the symmetry discrete along the length of the 
cylinder. 

When using a DMRG unit cell of 5m flue, the wavelength 

OiirB 

of the CDO is forced to be A m = 5m-^. As the circumfer¬ 
ence L is changed, A m deviates from the intrinsic wavelength 
Acdo of the CDO, and the stripes are forced to compress or 
stretch. Consequently we expect the energy density of the 
CDO depends parabolically on L, via E ~ (A m — Acdo) 2 , us 
found in Fig. |5(a)| For w = 0 we verified that this feature is 
found both for m = 3, L ~ 201b and the larger circumfer¬ 
ence m = 4, L ~ 27£b- The optimal Acdo can be determined 
from the minimum of the parabola, and depends on w at the 
level of 10 %, with larger w preferring smaller Acdo- 

As shown in Fig. |5(c)[ the CDO has remarkably close en¬ 
ergy (AE’/flux ~ 10 “^) to the RR 3 phase. In fact, for 
L rsj 20 £b and w = 0, the CDO is the true ground state. 
Increasing w favors the RR 3 state; our finite size numerics 
points to a transition around w ~ 1 £b • Our numerics on the 
infinite cylinder shows a sharp first order transition between 
the RR 3 and CDO phases. 

While Acdo ~ 4.75^^ is rather close to Ahf, the CDO 
we find is a highly entangled state rather different from the 
naive HF ansatz. In the HF ansatz for a stripe, the orbital 
occupation alternates between v = 1 and v = 0 , effec¬ 
tively forming decoupled stripes of IQH. Because the IQH 
has chiral edge modes, the correlations are algebraic along the 
length of the stripe. As shown in Fig. |5(b)[ we find the CDO 
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FIG. 5. (a) Energies for the RR3 states (red and blue triangles) as a function of L and CDO (green circles) as a function of wavelength 
A, at with zero well-width. For data points denoted by triangles, the energies are computed while enforcing translational invariance of the 
wavefunctions. The two ground states |^i/ T ) are metastable, and their energy approaches Err « —0.13717 in the thermodynamic limit. For 
data points denoted by circles, a charge density wave of wavelength A is imposed along the length of the cylinder (illustrated in the inset). 
Here the CDO minimum Ecbo dips below Err, which suggests that at w = 0, the quantum Hall phase is unstable to formation of CDO. (b) 
Orbital (dots) and real-space (line) density profile of the CDO minima at A ~ 4.8£b- Inset shows the Fourier transform of the orbital density, 
with an exponential decay, (c) Energy splitting (per flux) between the RR3 phase and the CDO phase (Ecbo — Err), as a function of w. The 
data shows that increased well-width w tends to stabilize the Read-Rezayi phase. 


orbital occupation instead forms a nearly perfect sine-wave, 
with higher Fourier components falling off exponentially in 
k. This implies that the Greens function decays exponentially 
along the length of the stripes, suggesting interaction effects 
have gapped out the chiral edge modes. This is consistent with 
theoretical predictions that the HF smectic order is unstable 
towards bubble (e.g. Wigner crystal) formation.® 61 We can 
observe the onset of a Wigner-crystal phase in the density- 
density correlations S( r) = (p(r)p(O)) shown in Fig. [6] 
Within the same stripe, the correlations are consistent with 
a crystal with one-electron per unit cell. Between neighboring 
stripes the correlations are much weaker, but suggest a body- 
centered rectangular structure. Thus while we cannot deter¬ 
mine if true long-range order develops between the stripes in 
the thermodynamic limit, the correlations are consistent with 
a body-centered rectangular lattice of electrons. 


B. Effect of Landau level mixing on the competition between 
the Read-Rezayi and re-entrant integer quantum Hall effect 

Since the ground states at v = ^ and are fully spin- 
polarized, in the absence of LL-mixing the two are related by 
particle-hole symmetry. To account for the discrepancy seen 
in experiment we now include LL-mixing, whose strength is 

parameterized by the ratio of Coulomb to cyclotron energies, 

2 

k = hjj c . We include the effect of LL-mixing by per¬ 
forming DMRG in the full space of spin-polarized N = 0,1, 2 
LLs . 34 Adding LL-mixing is computationally expensive; each 
data point takes several weeks. 

We first outline the logic of our approach. Since it appears 
that there is a first-order transition between the RR 3 states 
and CDO, the relevant task is not to compare the properties 


of RR3 vs. RR3 as LL-mixing is introduced, but rather to 
compare the energy of the RR3 state to the competing CDO 
at the same filling. In particular, near a first order transition the 
quasiparticle gap and correlation functions needn’t be related 
to the extensive energy splitting between the competing phases 
which actually drives the transition. Thus, we must compute 
the ground state energy density of the RR3 states and CDO 
states, at both filling = and as LL-mixing is added. 

As LL-mixing is increased the CDO wavelength Acdo 
changes slightly, so it is important to repeat the analysis of 
varying the circumference L. In Fig. [7] we compute the en¬ 
ergy difference between the CDO and RR 3 states at k = 0.4, 
w = 3 £3 for various circumferences. We find an asymmetry 
between the two fractions. At v = LL-mixing increases 
the energy of the CDO state relative to the RR 3 state, thus 
stabilizing the RR 3 quantum Hall plateau. For our data at 
v = however, the RR 3 -numerics become much more un¬ 
stable in the presence of LL-mixing and it is difficult to con¬ 
verge our numerics with the accuracy we require (with LL- 
mixing we have data up to x = 14000). Our uncertainty in the 
DMRG ground state energy is indicated by error bars, which 
we determined as follows. DMRG produces a variational en¬ 
ergy E(x) which depends on the DMRG bond-dimension x 
as E(x) = Eo + where Eq is the true ground state 

energy . 7 A fit to this form is done with and without the largest 
available bond-dimension, and we consider the change in the 
fitting parameter Eq to be the uncertainty. While arguably a 
flawed measure, by this definition the change in the CDO- 
RR 3 energy splitting is inconclusive at ^. However, our best 
data point is consistent with a scenario which the CDO-RR 3 
splitting changes significantly less than the CDO-RR 3 split¬ 
ting at 4 ?, if at all. 

The data in Fig. [7] is taken at LL-mixing n = 0.4, smaller 
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x (along cylinder length) 

FIG. 6. The real space density-density correlation function S(r) = 
(p(r)p(O)) of the CDO phase, at v — 44, for w = 1 £b, and 
L = 26.5 £b- The data is scaled by the mean density-squared 
( 27r ^ 2 ) 2 , and the point r = 0 is taken to lie at the center of the 
high-density stripe. The horizontal axis is along the length of the 
cylinder; the vertical axis is around the circumference of the cylin¬ 
der. At this circumference the unit cell along the cylinder contains 
4 • 5 flux, corresponding to eight electrons per stripe. The intra¬ 
stripe correlations have seven local maxima plus a correlation hole. 
The nearest-stripe correlations have weaker maxima which are out 
of phase with the intra-stripe maxima, which suggests there is a ten¬ 
dency towards forming a body-centered rectangular lattice of single 
electrons. 

than that of the experiments in Ref. [16| (k ~ 1). At least per- 
turbatively the change in splitting should scale linearly with 
ft, implying the true energy difference is even greater than re¬ 
ported here, though we suspect ft = 1 is beyond the perturba¬ 
tive regime. We only present data at ft = 0.4 because we are 
unable to converge the ft = 44 RR 3 state at larger ft, appar¬ 
ently because much larger bond dimension becomes required 
(the RR3 state at ft = 4= does not have as severe a depen¬ 
dence). Note that while the computational difficulty of RR3 
increases relative to RR3 there is not much relative difference 
in the observable properties like the correlation length, and no 
indication of a diverging correlation length. 


VI. SUMMARY AND DISCUSSION 

We numerically simulated the fractional quantum Hall ef¬ 
fect at a filling factor v = 44, 44 using the infinite density 
matrix renormalization group and exact diagonalization meth¬ 
ods. Our simulations include realistic Coulomb interaction 
appropriate for GaAs quantum wells of finite width. In the 
absence of Landau level mixing, the topological properties 
of the ft = 44 ground state are consistent with the k = 3 
Read-Rezayi phase. The ground state remains spin-polarized 
for a range of well-widths, even as the Zeeman splitting van¬ 
ishes. The lowest energy (spin-polarized) charged excitation 
was identified with a non-Abelian Fibonacci anyon, which 
supports universal braid statistics. 


FIG. 7. The effect of Landau level mixing on energy splitting (per 
flux) between the CDO and the Read-Rezayi states. The data is given 
for w = 3 £b and ft = 0.4, where ft is the ratio of Coulomb to cy¬ 
clotron energy. For reference, the dashed line denotes the energy 
splitting in the absence of LL-mixing. For v — 44, the splitting in¬ 
creases with the inclusion of LL-mixing, further favoring the RR3 
phase. For v — 44 the data is more ambiguous due to our numerical 
uncertainty in the energy of the RR3 phase, which becomes signifi¬ 
cantly less stable with the addition of LL-mixing. We conclude that 
LL-mixing stabilizes the RR3 phase at ft = 44. 


Our findings reveal several properties which could help 
make the 12/5-plateau suitable for a bulk implementa¬ 
tion of universal measurement-only topological quantum 
computing.® First, our numerically calculated gap of 2 K sug¬ 
gests that the 12/5 Fibonacci state is not intrinsically much 
more delicate than the 5/2 state. The smaller 80 mK ob¬ 
served in experiment may be related to slowly varying disor¬ 
der potentials, and hence there is some hope this gap could 
be engineered. Second, of the two types of charge ±e/5 
quasiparticles, the Fibonacci ±e/5 r particle has the lower en¬ 
ergy. Furthermore, the ±e/5r particles and ±e/5 particles 
have a qualitatively different charge profile; the Fibonaccis 
are much more tightly localized at their center (see Fig. [3}. 
This implies that an electrostatic pinning potential will sta¬ 
bly trap a Fibonacci anyon; in contrast, if it was the Abelian 
=be/5 which had a much lower energy when pinned, a trapped 
=be/5r would spontaneously eject the neutral r, which would 
fly off in an uncontrolled manner. This also provides a con¬ 
venient way to initialize a lattice of Fibonaccis with a known 
fusion tree: adiabatically turning on two nearby pinning po¬ 
tentials +V, —V would usually generate a e/5r, — e/5r pair 
that fuse to the identity, which can then be separated. In ad¬ 
dition, rather than using an edge-interferometer to detect fu¬ 
sion outcomes, one can try a more naive approach: bring two 
charged Fibonaccis together and “look” at them. Since we 
find the charge distribution depends qualitatively on the fu¬ 
sion outcome, measuring, for instance, the local quadrupole 
moment or the magnetization density detects the fusion out¬ 
come. Thus, given sufficient control over local pinning po¬ 
tentials, we have the basic required ingredients: 1) charged 
Fibonaccis can be stably trapped by electrostatic pinning po¬ 
tentials 2) an array of Fibonaccis can be pairwise created with 
a known fusion outcome 3) the charged Fibonacci and Abelian 
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particles have qualitatively different charge profiles, so fusion 
outcomes can be detected by any probe sensitive to the local 
charge or magnetization distribution. 

The full spin polarization of the ground state and the large 
charge gap we obtain (0.017 « 2K) are encouraging but 
not in complete agreement with experiments The esti¬ 
mated gap in Ref. [16] is « 80 mK, and Ref. [49] detected a 
spin transition upon tilting the field. While it may be possible 
that the lowest charge excitation is actually a skyrmion and/or 
strongly renormalized by the LL-mixing, a better theoretical 
understanding of activated transport and spin polarization in 
the exotic N = 1 LL plateaus is also desired. As discussed in 
Ref.[62]in the context of v = 5/2 plateau, the size of the quasi¬ 
particles (d ~ 1 Mb) is comparable to the length scale of dis¬ 
order arising from the remote ionized donors. In this regime 
there is a larger tunneling amplitude across saddle points of 
the disorder potential; it would be interesting if these ampli¬ 
tudes could be numerically estimated using the single-anyon 
DMRG. 

Perhaps most intriguingly, we find an exceptionally close 
competition between CDO and the RR 3 phase, the former be¬ 
ing the likely origin of the RIQH phase experimentally ob¬ 
served at v = Our numerics show that with increased 
width w the RR 3 phase is preferred over the CDO; it is ad¬ 
vantageous to fabricate the quantum well as wide as possible 
to stabilize the RR 3 phase. At the same time, with increased w 
the 1 st excited subband LL also comes down in energy, cross¬ 
ing with the N = 1 LL at w ~ 3.8 £b> which puts an upper 
limit on w. Finally, Landau level mixing increases the energy 
of the CDO state relative to the RR 3 at v = ^, which would 
explain the observed asymmetry between the v = ^ and ^ 
plateaus, giving further confidence to the numerics. 

Note: During the final preparation of this work, we learned 
of overlapping results by W. Zhu et al.p^and a study of Lan¬ 
dau level mixing on the sphere by Pakrouski et al.® 
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Appendix A: Further characterization of the Read-Rezayi phase 
via entanglement 


We discuss in details the various entanglement measures to 
identify the RR 3 phase at filling v = ^. The data presented 
here are are computed assuming full spin-polarization and no 
Landau level mixing, i.e., a single N — 1 Landau level at 
| filling. (The data at w = 0 are taken while enforcing a 
uniform density to stabilize the RR 3 phase. It appears rather 
“noisy”, possibly due to a CDO instability.) 

We first return to the momentum polarization for the ground 
states \Clt) and \Cl T ). As alluded to in the Sec. |n| the mo¬ 
mentum polarization reveals the shift, chiral central charge c, 
and topological spin h a . We numerically calculate the Berry 
phase Ut-o, of performing a 2i r twist on the left half of the 
cylinder for ground state |fi a ): comparing to the theoretical 
predictiorP^ 


U T -a = exp 


2m 



S _ _ Sh_ l 2 \ 
24 2irh ) 



(Al) 


The ellipsis denotes term exponentially suppressed with cir¬ 
cumference L , and r\ H is the “Hall viscosity”,^ related to the 
shift via r] H = | ^(For the RR 3 phase, the shift is com¬ 


puted at v = | in the 1 st LL.) The formula given in the main 
text is the “logarithm” of Eq. ( |A1| ), and in principle ambiguous 
modulo 1. In practice, because Ur-a is a continuous function 
of the circumference, one can easily extract the S by taking 
the slope and compute the combination h — ^ modulo 1 . 

Figure 8 (a) shows the momentum polarization for L = 
15-22^5, using various range of thicknesses w. From the 
slope of the data, we extract S = 0 db 0.1, in excellent agree¬ 
ment with the RR 3 phase at zero shift. Assuming the shift is 
exactly zero, the data shows the residue h a — ^. These values 


corroborate well with those of the RR 3 phase with h\ = 0, 
h T = — | and chiral central charge c = 1 — | = —|, shown 
as the pair of dashed lines. 

Finally, we also study the entanglement entropy of the 
ground states, given by S = Tr[— pL logpzj. The entangle¬ 
ment entropy is expected to scale with the circumference as 
S a « sL — 7 a , where 7 a is a constant called the “topologi¬ 
cal entanglement entropy” (TEE) associated with any on type 
a> 37]66| Th e TEE is given by 7 a = log V — log d a , where d a are 
the quantum dimensions of the any on type a, and V is the to¬ 
tal quantum dimension of the system given by V 2 = s ^ Ja d 2 a . 
Figure [ 8 ] shows the entanglement entropy as a function of L. 
Unfortunately the data S t / r suffer from very strong finite size 
and finite entanglement effects, it is not possible to extract 
7i/t with any meaningful degree of certainty (the slope s is 
sensitive to microscopic details and non-universal). However, 
the difference S T — St = 71 — 7 r is universal, predicted to 
be log((/?) where ip = 1+ 2 V ^ is the golden ratio. While the en¬ 
tropy data alone cannot confirm the existence of a Fibonacci 
anyon, it is nevertheless consistent with the RR3 phase. 



(b) 6.0 



0.0 — 1 — 1 — 1 —*—1—*—1—*—1—'—1—'—1—■—1— 
15 16 17 18 19 20 21 22 

lub 


FIG. 8. (a) Momentum polarization for the ground states \Qt/ T ). The 
dashed lines are the predicted values of h a — ^ for the two ground 
states, (b) Entanglement entropies St and *SV, along with their dif¬ 
ference S T — St- The dashed line at log(c^) « 0.48 is predicted for 
the RR3 phase. 



FIG. 9. Spin correlator (Sq S k ) for L — 2 1 £b, w = 3 £b at filling 
v t/4- — 1. Solid curves shows the data at various bond dimensions 
X, converging to the expected correlator given as the dashed curve. 
This is a signature of spontaneously broken spin-rotational symme¬ 
try, indicative of a fully spin-polarized ground state at v — ^. 
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Appendix B: Stability of the Read-Rezayi phase against 
short-range perturbations 


The ground state at v = 2/5 in the TV = 0 LL is the spin- 
polarized Abelian hierarchy (AH) or composite fermion state, 
and an obvious competitor at v = 12/5. Here we examine 
the stability of the RR 3 phase at v = 12/5 as we perturb the 
TV = 1 LL projected Coulomb interaction by the short-range 
V\ pseudopotential. 

As the node of TV = 1 LL wavefunctions softens the inter¬ 
action, we expect that adding AVi > 0 will drive the RR3 
phase back into the Hierarchy phase. This is indeed what hap¬ 
pens, Fig. 10 In terms of the ratio V1/V3, the transition be¬ 
tween the RR and hierarchy phase is about 1/3 of the way 
between the Coulomb N m l and TV = 0 LLs. AVi is a 
‘best-case’ perturbation for the Hierarchy state, as a typical 
real-space potential will be distributed across all V m . In fact, 
we have verified that projecting a V(r) = V 2 £(r) interaction 
into the TV = 1 LL favors a CDO phase, not the Hierarchy 
state. Linearly extrapolating the energy of the Hierarchy state 
to the Coulomb point, we obtain a splitting between the RR3 
and Hierarchy state of about A E ^ 1.5 x 10 -4 per flux at 
well-width w = 2 £b • 

We directly probe the stability of the entire low-energy 
spectrum of the RR 3 phase upon varying V\ using exact di- 
agonalization (ED) [Fig. ED- We assume complete spin po¬ 
larization and no Landau level mixing. The energy spectrum 
of the system is resolved as a function of pseudomomentum 
K El and the RR 3 phase is characterized by a two-fold degen¬ 
erate ground state in K = 0 sector. Additionally, as discussed 
in Sec. [IIJ there are five copies of those that are related by the 
center of mass translation and can be factored out. 

In Fig. [IT] we show the phase diagram of the system as V\ 
is modified, for zero width (left) and w/£b = 4 (right). We 
compute 10 lowest energies per momentum sector of the sys- 



FIG. 10. Transition between the RR 3 and Hierarchy states as the 
Haldane pseudopotential AVi is added to the Coulomb interaction. 
In the cylinder geometry the transition is first order. For perspective, 
in the TV = 0 LL, V 1 /V 3 = 1.6; in the TV = 1 LL, V 1 /V 3 = 1.3; 
at the observed transition V 1 /V 3 — 1.43. Linearly extrapolating the 
energy of Hierarchy state to the Coulomb point, we obtain a splitting 
of Eau — Err % 1.5 x 10 _4 /flux. Data is taken from DMRG at 
L = 20fe, w = 2fe. 



FIG. 11. Phase diagram as a function of modified Vi pseudopoten¬ 
tial, for zero width (left) and w/£b = 4 (right). Data is obtained by 
exact diagonalization of the system with 25 flux quanta through the 
hexagonal unit cell. We compute the lowest 10 energies per momen¬ 
tum sector and plot them relative to the average energy of the system 
E — E — E awg . Black symbols denote levels in K = 0 sector. For 
zero width (left), we identify the Z3 ground state degeneracy in a nar¬ 
row shaded region around the Coulomb point SV 1 = 0. In this region 
the ground state has large overlap O with the model Z3 wavefunc- 
tion (inset). At width w/£b — 4 (right), the Z3 phase widens and 
becomes more robust. It is surrounded by the hierarchy/composite 
fermion state for larger positive 5V 1 , and several charge density or¬ 
dered phases for negative SV 1 . 


tern with 25 flux quanta through the hexagonal unit cell. The 
energies are given in units of E c , and for clarity we plot them 
relative to the average energy of the system E = E — E awg . 
Black symbols denote the levels belonging to K = 0 sector. 
For zero width (left), we identify the RR 3 ground state de¬ 
generacy in a narrow shaded region around the Coulomb point 
SVi = 0. In this region the ground state also has large overlap 
O with the model RR 3 wavefunction (inset). Because of the 
strong mixing of four lowest energy levels with K = 0 around 
the Coulomb point, we define the overlap O as a sum of sin¬ 
gular values of the 4 x 2 overlap matrix Oij = (V4xactl^z 3 )’ 
i = 1,..., 4, j = 1 , 2. At larger width w/£b = 4, the RR 3 
phase widens and becomes more robust as the two quaside¬ 
generate K = 0 levels become better separated from the rest 
of the spectrum. By inspecting the level degeneracy and com¬ 
puting the overlaps, we also deduce that the RR 3 phase is sur¬ 
rounded by the hierarchy/composite fermion state for larger 
positive SVi, and several charge density ordered phases for 
negative SV The phase CDO I was identified with a stripe 
in Ref. |27| The estimate of critical V\ for the transition into 
the Hierarchy state is in agreement with DMRG estimate in 
Fig- [TO] 



























11 



FIG. 12. Signature of spontaneous spin rotation symmetry breaking. (Left) The entanglement spectrum of the spin-full system at 

sorted by charge sector (Q^, Q^). (Right) The entanglement spectrum for a spinless v “ | system for comparison. Both spectra are taken at 


L = 2U b ,w = 3£ B - 


Appendix C: Spin polarization 


We first clarify the expected signatures of spin-polarization when explicitly preserving S z invariance at filling iA = = 

If I^RRs ^presents a RR 3 state spontaneously polarized in the XY plane at angle 6 , we force the state into the superposition 


l*> = / d 

Jo 


d® I^)rr 3 


(Cl) 


This state has infinite bipartite entanglement, so cannot be represented exactly by a matrix product state (MPS). However, we 
expect that the (S + {r)S~ (r')) correlations will be large and nearly constant out to a distance which will depend on the MPS 
bond-dimension x, after which it will decay exponentially. (We define 5+ = cjc^, and S~ to be its Hermitian conjugate.) This 
effect (in orbital space) is shown in Fig. [9] for L = 21£b, w = 3£b- In the x °° limit, we expect the correlation function to 
approach the number-number correlator ^ (n 0 n/c) for the fully spin-polarized RR 3 phase (shown in dashed line). 

The entanglement spectrum provides further evidence. Recall that the entanglement spectrum can be sorted by the U(l) 
charges (Q^, Q^). In Fig. 


12 


we show the entanglement spectrum for a 5 x 5 grid of (Q^, Q^) sectors. For constant Q, the line 
= Q contains many copies of the charge Q RR 3 spectrum, up to some cutoff \Q^ — Q^\ < *S™ ax (x) that depends on x- 
This is because in order to support long-range XY correlations, there must be very large fluctuations in S z across any bipartition 
of the system. The entanglement spectrum of a symmetry broken phase was discussed in Ref. [ 68 ] where it was shown that the 
spectrum contain a “tower of states” associated with the broken symmetry; the many identical copies of RR 3 we find is this 
tower of states. 
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FIG. 13. Real-space structure factor S( r) = (p(r)p(O)) / (p) 2 of the v — 2? Read-Rezayi and CDO states at L — 20 £b, w = 11b . The 
degenerate RR3 states are mostly identical which is a necessary condition for topological order. The RR3 states are also nearly rotationally 
symmetric in contrast to the striped structure of the competing CDO phase; the residual anisotropy is presumably due to the finite circumference 
of the cylinder. 


We also note that similar computation was perform for v 


^ and we also found evidence for spontaneous spin-polarization. 


Appendix D: Structure factor of the Read-Rezayi phase 


In Fig. 13 we show the structure factor S(r) for the RR 3 states and the CDO state. The ground states \Qj) and |O r ) have 


nearly identical structure factor, a requisite for topological order in the RR 3 phase. At L = 20 £b and w = 1£b, the energies 
of the RR 3 and CDO states are within 2 x 10 -5 of one another. However, despite having such close energies, the RR 3 states 
and CDO state have much different structure factors. Furthermore, while the CDO state has higher energy than the RR 3 states, 
it remains metastable in our DMRG simulations indicative of a 1 st order transition. 


r 



FIG. 14. The effect of finite width w/£b on the hexagonal toms 
threaded by = 25 flux quanta (left) and N$ = 30 (right). Energy 
spectrum is obtained by exact diagonalization and plotted relative 
to the average energy of the system. Black symbols denote levels 
belonging to K = 0 sector. 


Appendix E: Exact diagonalization: effects of finite well-width 
and torus aspect ratio 


We showed in Sec. [IJ and [V] that finite width stabilizes the 
homogeneous fluid phase. In Fig. [l4jleft) we systematically 
study the effect of well width on the energy spectrum (in par¬ 
ticular, the topological ground state degeneracy) via exact di¬ 


agonalization. We consider a hexagonal torus threaded by 
N<$> = 25 and N<$> = 30 flux quanta. 


For a smaller system (TV$ = 25), the topological degener¬ 
acy is not well-resolved for zero width due to the mixing with 
a higher level belonging to K = 0 sector (black symbols). In 
this case, the main effect of non-zero w/£b is to lift the spuri¬ 
ous level in energy, leaving a robust two-fold degenerate man¬ 
ifold of ground states. The gap separating these two quaside¬ 
generate states from the rest of the spectrum further widens 
as w/£b is increased. At the same time, the overlap with the 
model Read-Rezayi state slightly increases as a function of 
width (not shown). For a larger system (7V$ = 30), the two¬ 
fold degeneracy appears to be present already for zero width, 
but it gets better resolved for moderate widths w/£b ~ 1.7. 


In Fig. [l5jtop) we show the effect of changing the geome¬ 
try of the torus unit cell. We consider a rectangular L x x L y 
torus in this case, whose area is fixed due by the condition 
L x L y = 27r£gN$. By changing one of the linear dimen¬ 
sions of the torus (L x ), we can drive a transition between the 
homogenous phase and the CDO. In Fig. F5 L x « Ills cor¬ 
responds to an isotropic torus (L x = L y ) where the ground 
state is approximately two-fold degenerate and belongs to the 
Read-Rezayi phase. Beyond L x « 19^, the system evolves 
towards another, much deeper, energy minimum, which was 
identified in Sec. |V| with the CDO phase. The difference in 
ordering between the two phases is captured by the pair corre¬ 
lation function g(r) = N L ^ v _^ (5(r — + Rj)) shown in 
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FIG. 15. The effect of changing the geometry (one linear dimension 
L x of the toms) on the energy spectrum (top) and the ground state 
(bottom row) of the system. Energy spectrum is obtained by exact 
diagonalization of the rectangular toms with N$ — 20 flux quanta, 
and black symbols denote levels belonging to K = 0 sector. Varying 
L x induces charge density order in the ground state, as seen in the 
pair correlation function g(r) for several values of L x (bottom). 


Fig.fHJbottom). Unfortunately, because the area of the torus 
must be preserved as we change L x , this implies that corre¬ 
lations along the ^/-direction in the ground state of the sys¬ 
tem at L x & 19 £b are artificially truncated because of small 
L y . Therefore, the CDO ground state in this case is likely not 
faithfully reproduced due to finite-size effects and has signif¬ 
icantly less entanglement than what we found by DMRG in 
Sec. [VI 


Appendix F: Defect-DMRG 

Here we briefly discuss the method used to generate pinned 
anyons. The iMPS ansatz for a single any on excitation was 


discussed in Ref. 0Tj Because the anyons are charged, in the 
absence of a pinning potential the single-anyon states form a 
Landau level. It is convenient to study the anyons in the Lan¬ 
dau gauge so that anyon a forms a plane wave with momen- 

i 1 2 3 4 

turn k around the cylinder and is localized near x oc k along 
the length. This choice allows us to conserve momentum in 
the DMRG simulations (interestingly, k is actually fractional 
because of the topological spin of the anyon). After DMRG 
we obtain single-anyon states | k) a \ just like a Landau level, 
their energy is independent of k , and different | k) are related 
by a magnetic translation T x along the length of the cylinder. 
The energy E a reported in the main text is the energy of \k) a 
relative to the ground state, combined with the interaction of 
a point-charge with a neutralizing background. 

We caution that we have not performed finite size scaling 
of the gaps. If one is familiar with gap calculations on the 
sphere, this might seem problematic: on a sphere with N elec¬ 
trons, the corrections are of order 1/N and it is imperative to 
extrapolate gaps in 1/N. However, the absence of curvature 
on the cylinder leads to a much more favorable scaling. One 
can show that if the quasiparticle density is bounded by an 
exponential tail, our estimate of the gap will converge expo¬ 
nentially quickly in the cylinder circumference to its Coulomb 
value. (Strictly speaking, for a Coulomb interaction the quasi¬ 
particles will have 1 /r 5 6 tails, but this will lead to weak cor¬ 
rections.) As a benchmark, if we apply our procedure to cal¬ 
culate the energy of a quasihole in the integer quantum Hall 
effect at the same circumference (L = 21 £b) and precision as 
the results quoted in the text, we obtain a quasihole energy of 
E « 1.2531; the exact result is ~ 1.2533. 

To pin the anyons in 2D we need to introduce a pinning 
potential c/)(x,y), but it is too expensive to work in the full 
Hilbert space without conserving momentum. In the limit of 
a weak pinning potential we can use first-order degenerate 
perturbation theory and project </> into the variational space 
spanned by the anyonic Landau level {|fc) a }. We compute the 
many-body matrix elements = (k'\ \ k) a using stan¬ 

dard MPS techniques, diagonalize iT eff , and compute the real- 
space density of the lowest energy state |fc) a . 

The total energy of the pinned excitation is the energy E a 
reported in the main text plus the potential energy H Qff from 
the pin. For a width 2 Ib Gaussian, the pin lowers the en¬ 
ergy of the non-Abelian |r excitation by 85% more than the 
Abelian for a width 6 is Gaussian, the difference is 25%. 
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